Critical behavior of the conductivity of Si:P at the metal-insulator transition under 

uniaxial stress 
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We report new measurements of the electrical conductivity a of the canonical three-dimensional 
metal-insulator system Si:P under uniaxial stress S. The zero-temperature extrapolation of a(S, T — + 
0) ~| S — S c | M shows an unprecidentedly sharp onset of finite conductivity at S c with an exponent 
H = 1. The value of differs significantly from that of earlier stress-tuning results. Our data show 
dynamical a(S, T) scaling on both metallic and insulating sides , viz. a(S, T) = cr c (T) ■ T{\ S — S c \ 
/T y ) where o c {T) is the conductivity at the critical stress S c . We find y — 1/zu = 0.34 where v is 
the correlation- length exponent and z the dynamic critical exponent. 
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Quantum phase transitions have become of steadily 
increasing interest in recent years (lj. These continuous 
transitions ideally occur at temperature T — where 
quantum fluctuations play the role corresponding to ther- 
mal fluctuations in classical phase transitions. In partic- 
ular, certain types of metal-insulator transitions (MIT) 
such as localization transitions have been studied exten- 
sively. Experimentally, the MIT may be driven by an 
external parameter t such as carrier concentration N, 
uniaxial stress S, or electric or magnetic fields. Gener- 
ally, electron localization might arise from disorder (An- 
derson transition) or from electron-electron (e-e) inter- 
actions (Mott-Hubbard transition) 0. In Nature, these 
two features go hand in hand. For instance, the disorder- 
induced MIT occurring as a function of doping in three- 
dimensional (d = 3) semiconductors where the disorder 
stems from the statistical distribution of dopant atoms in 
the crystalline host, bears signatures of e-e interactions as 
evidenced from the transport properties in both metallic 
H and insulating regimes pj. This makes a theoretical 
treatment of the critical behavior of a MIT exceedingly 
difficult. Even for purely disorder-induced transitions, 
the critical behavior of the zero-temperature dc conduc- 
tivity, cr(0) ~| t — t c | M where t c is the critical value of 
t, is not well understood. Theoretically, fi is usually in- 
ferred from the correlation-length critical exponent v via 
Wegner scaling [i — v(d—2). Numerical values of v range 
between 1.3 and 1.6 |g,D. 

Experimentally, it has long been suggested that the 
critical behavior of the conductivity falls into two classes: 
/i w 0.5 for uncompensated semiconductors and /i sa 1 for 
compensated semiconductors and amorphous metals [Q. 
However, there appears to be no clear physical distinc- 
tion between these materials that would justify different 
universality classes. While many different materials were 
reported to show /i w 1, the exponent /i w 0.5 was largely 
based on the very elegant experiments by Paalanen and 
coworkers iBHlQl i where uniaxial stress was used to drive 



an initially insulating uncompensated Si:P sample metal- 
lic. This allows to fine-tune the MIT since the stress can 
be changed continuously at low T thus eliminating geom- 
etry errors incurring when different samples are employed 
in concentration tuning the MIT. 

As always when dealing with critical phenomena, the 
range of critical behavior is a source of controversy. A few 
years ago we suggested JO] to limit the critical concentra- 
tion region in doped semiconductors on the metallic side 
of the MIT to samples where cr(T) actually decreases with 
decreasing T, i.e. the sample becomes less conducting 
when approaching the MIT. In doped semiconductors, 
cr(T) is nearly independent of T at the crossover concen- 
tration N cr , with a value a cr of a few times 10 ft cm , 
e.g. a cr w 40£! _1 cm -1 in Si:P. u(T) exhibits a negative 
temperature coefficient above N cr which is explained in 
terms of e-e interactions . Typically the critical region 
N c < N < N cr is within 10% or less of the critical con- 
centration N c . This eliminates a large number of studies 
purporting to show fi = 0.5 where actually only a few 
samples in the critical regime were investigated. Even 
the recent study on transmutation-doped Ge:Ga where 
fi = 0.5 was suggested, presents only three metallic sam- 
ples in the critical region below a cr ~ 10 51 _1 cm _1 [ [l3|| . 
An earlier study of a large number of Si:P samples showed 
that the conductivity exponent fj, changed from /i = 0.64 
for N > N cr w 1.1 N c to 1.3 for N c < N < N cr [0. 
On the other hand, sample inhomogeneities might affect 
the behavior very close to N c . For this reason, data for 
stress-tuning with stress close to the critical value were 
discarded in the earlier study, leading to /J, = 0.5 p|,|9|,|l2t . 
It is therefore absolutely necessary to perform additional 
stress-tuning experiments on Si:P with finely tuned stress 
values including data on the insulating side to check for 
the critical behavior. 

The notion of a quantum phase transition allows a sec- 
ond important aspect to be addressed, namely the inter- 
dependence of static and dynamic behavior. The dy- 
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namics is reflected in the finite-temperature behavior of 
critical quantities. Concerning the MIT in heavily-doped 
semiconductors, this point has not received much atten- 
tion from the experimental side. A first attempt was 
made fbl using the scaling function H 



a{t,T) = {t-t c yF{T/{t-t c f v ) 



(1) 



where z is the dynamic critical exponent. This relation 
is often referred to as dynamic scaling. Approximate dy- 
namic scaling was observed for Si:P on the metallic side 
of the MIT with t = N, yielding fi = 1.3 and z = 2.4 
|l4|| . On the other hand, the stress-tuning data ||] did 
not obey scaling Very recently, Bogdanovich et al. 
flBf demonstrated that conductivity data for Si:B un- 
der uniaxial stress obey very nicely the dynamic scaling 
on both metallic and insulating sides, yielding ji = 1.6 
and z = 2, while concentration tuning of cr(0) on the 
same system had suggested /i = 0.63 Jig] , This large 
difference is not understood at present. In this situa- 
tion, an examination of possible dynamic scaling of the 
canonical metal-insulator system Si:P appears of utmost 
importance in order to resolve the question of critical 
behavior and to appraise the possibly strongly different 
roles of stress and concentration in tuning the MIT. 

In this paper, we report on stress tuning of the MIT 
of Si:P by measuring the electrical conductivity down 
to 15 mK. We find by extrapolating to T — an un- 
precidently sharp onset of a(t, 0) which allows to unam- 
biguously extract |i « 1. In addition, dynamic scaling 
yielding z ~ 3 is found. The value of /z is in reasonable 
agreement with that derived from concentration tuning. 
We further demonstrate that stress tuning and concen- 
tration tuning lead to very different T dependencies of 
a. 

The samples were taken from the same Si:P crystals 
which have been employed previously Jll]]. Here we re- 
port on investigations on two crystals with N = 3.21 
and 3.43 -10 18 cm _3 , just below the critical concentration 
N c = 3.52- 10 18 cm -3 as determined |ll]] for our samples. 
Similarly grown samples with an even higher concentra- 
tion (N w 7 • 10 19 cm~ 3 ) showed no sign of P clustering 
as investigated with scanning tunneling microscopy [fl7| . 
The samples were cut to a size of ~ 15x0.8x0.9 mm 3 and 
contacted with four Au leads by spark welding, with the 
voltage leads ~ 6 mm apart. The sample was mounted 
in a 4 He-activated uniaxial pressure cell equipped with a 
piezoelectric force sensor. The stress was applied along 
the [100] direction which was the most elongated dimen- 
sion of the sample. The stress was determined from the 
ratio of the area of the cell base plate and the sample 
cross section. Calibration of the cell showed a linear in- 
crease of force with pressure applied at room tempera- 
ture to gazeous He, with no hysteresis between increasing 
and decreasing pressure. The cell, incorporating a ther- 
mal shield, was tightly screwed to the mixing chamber 



of a dilution refrigerator. During one run a thermome- 
ter was attached to the sample showing that temperature 
deviations to the main thermometer directly mounted at 
the mixing chamber were less than 0.5 mK at the lowest 
measuring temperature of 15 mK. The conductance was 
measured with a LR 700 resistance bridge at 16 Hz. 
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FIG. 1. Conductivity a of a Si:P sample with P concen- 
tration N = 3.21 • 10 18 cm~ 3 versus y/T for several values of 
uniaxial stress. From top to bottom: S = 3.05, 2.78, 2.57, 
2.34, 2.17, 2.00, 1.94, 1.87, 1.82, 1.77, 1.72, 1.66, 1.61, 1.56, 
1.50, 1.41, 1.33, 1.26, 1.18, 1.00 kbar. Solid lines are connect- 
ing the very finely spaced individual data points. 

Fig. 1 shows the electrical conductivity <r(T) of sam- 
ple 1 (N — 3.21-10 18 cm -3 ) for uniaxial pressures between 
1 and 3.05 kbar. The data are plotted vs. \JT which is 
the T dependence expected due to e-e interactions and 
indeed observed well above the MIT, a(T) = erg + rriy/T 
with m < ||. The smooth curves are in fact polygons 
connecting adjacent data points (see Fig. 2a for a set of 
actual data points) . Under uniaxial stress between 1 and 
2.57 kbar the er(T) curves evolve smoothly from insulat- 
ing to metallic behavior with m > 0, and cr(T) becomes 
nearly independent of T with a value a cr ~ 12 f2 cm -1 
at ~ 2.7 kbar. For larger stress <r(T) passes over a shallow 
maximum signaling the crossover to m < 0, as observed 
with concentration tuning p8[ . It is interesting to note 
that a cr (S) rs 0.3a cr (N), thus severely limiting the crit- 
ical region. Our data do not exhibit the precipitous drop 
of cr(T) below ~ 40 mK for pressures closest to the MIT, 
in distinction to the earlier stress-tuning work on Si:P 
extending to 3mK [^]||. Instead, our cr(T) data exhibit 
a T dependence that varies only gently with stress. 

Closer inspection shows that the data near the MIT 
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FIG. 2. (a) Low-temperature data of a of Fig. 1 in the 
immediate vicinity of the metal-insulator transition plotted 
against T 1 ' 3 . Dashed line indicates the conductivity at the 
critical stress (see text), (b) Extrapolated conductivity <r(0) 
for T — > versus uniaxial stress S for two P concentrations 
N = 3.21 and 3.43 • 10 18 cm -3 (closed and open circles, re- 
spectively). The inset shows earlier <r(0) versus S — S c data 
(triangles) from Ref. [8] in comparison to our data for sample 
1 (circles). 



are actually better described by a T 1 / 3 dependence for 
low T as can be seen from Fig. 2a for a few selected pres- 
sures in the immediate vicinity of the MIT. cr(0) obtained 
from the T 1 / 3 extrapolation to T = is shown in Fig. 2b, 
together with data for sample 2 (N = 3.43 • 10 18 cm~ 3 ) 
(see Fig. 3 for a(T) of this sample for a few represen- 
tative uniaxial pressures). c(0) is plotted linearly vs. 
S, yielding S c — 1.75 kbar for sample 1 and 1.54 kbar 
for sample 2. Note that the critical stress S c is quite 
well defined, as <r(0) breaks away roughly linearly from 
zero within less than 0.1 kbar. Applying our criterion for 
the critical region, the analysis should be limited to data 
with a < u C r ~ 12 fi^cm -1 . In this range the critical 
exponent /i is 0.96 and 1.09 for sample 1 and 2, respec- 
tively, fi ~ 1 are found also when the more conventional 
^JT extrapolation is employed as can be inferred from 
Fig. 1. This behavior contrasts with the earlier stress- 
tuning data p reproduced in the inset of Fig. 2b, where 
appreciable rounding close to N c is visible as compared 
to our samples when plotted against S—S c (see also Q). 
However, those <r(0) data between 4 and 16 f2 -1 cm~Hire 
compatible with linear dependence on uniaxial stress. 
Fig. 3 shows a(T) of sample 2 for a range of selected 
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FIG. 3. Comparison of the concentration dependence 
of tx(T) of Si:P (open symbols, from top to bottom: 
N = 3.60, 3.56, 3.50, 3.38- 10 18 cm- 3 , from Ref. [11]) and stress 
dependence of a(T) (closed symbols, N — 3.43 • 10 18 cm -3 , 
from top to bottom: S = 2.48, 2.11, 1.72, 1.32 kbar). 



uniaxial pressures, again applied along [100]. The over- 
all behavior is very similar to that of sample 1. The fact 
that a cr is the same for both samples is nevertheless sur- 
prising given the difference in S c . It has been suggested 
that tuning with S or N should yield the same critical 
exponents p|- p0Pl9| ]. The decrease of N c with uniaxial 
stress is attributed to the admixture of the more extended 
ls(E) and ls(T 2 ) excited states to the ls(Ai) ground- 
state of the valley-orbit split sixfold donor Is multiplet 
p9[ . However, comparison of er(T) for various S and N 
(Fig. 3) reveals that stress and concentration tuning lead 
to strikingly different T dependences of the conductiv- 
ity in the vicinity of the MIT. As the exact origin of the 
<t(T) behavior close to the MIT is unknown, we cannot 
offer an explanation for the different behavior which, of 
course, must arise from the change of donor wave func- 
tions under uniaxial stress. In this respect, experiments 
on similar samples for stress applied along different di- 
rections leading to different types of mixing among the 
states of the Is multiplet will be helpful. The fact that 
stress was applied to different directions in the previous 
and present study, i.e. [123] and [100], respectively, may 
well be one reason for the different behavior of a(T). 

We finally turn to the scaling behavior of a at finite 
temperatures using the data of sample 1. We employ the 
scaling relation [113] 
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a(t,T)=a c (T)F((t-t c )/Ty) 



(2) 



where a c (T) = a(t c ,T) is the conductivity at the critical 
value t c of the parameter t driving the MIT. This scaling 
relation is equivalent to Eq. (1), both are derived from 
the general scaling relation 



tr(t, T) = b-^-^T ((t - t c )b x /\ b z T) 



(3) 



where b is a scaling parameter. If the leading term to 
cr c (T) is proportional to T x , one obtains x — fi/vz and 
y = X/vz from a scaling plot. Fig. 1 and 2a show that 
a for 5 close to S c does not exhibit a simple power-law 
T dependence over the whole T range investigated. We 
therefore determine <J C (T) self-consistently in the follow- 
ing manner. The critical stress for sample 1 is taken from 
the above analysis as S c — 1.75 kbar. In order to ob- 
tain a c (T), we interpolate linearly between the two cr(T) 
curves for S = 1.72 and 1.77kbar. The resultant & C (T) is 
then fitted by the function <r c (T) = aT x (l + dT w ) with 
a = 6.01f2 _1 cm- x ,a; = 0.34, d = -0.202, w = 0.863, 
and T is expressed in K. Here the dT w term presents a 
correction to the critical dynamics. This & C (T) curve is 
shown in a dashed line in Fig. 2a. All a(S, T) curves with 
1.00 kbar < S < 2.34 kbar up to 800 mK are then used for 
the scaling analysis according to Eq.(2). The same proce- 
dure was repeated for other choices of cr c (T) between the 
two measured <r(T) curves embracing the critical stress 
with clearly less satisfactory results. 

Fig. 4 shows the resulting scaling plot of a(S, T)/a c (T) 
vs. | S — S c | / S c T y . The data are seen to collapse on a 
single branch each for the metallic and insulating side, 
respectively. The best scaling, as shown, is achieved 
for y = 1/zv = 0.34. Together with /u = 1.0 as ob- 
tained from Fig. 1 and assuming Wegner scaling v = fj, 
for d = 3, we find z — 2.94, which is indeed consistent 
with a c ~ T x l z ~ T 1 / 3 for T -> (see Fig. 2a). Alter- 
natively, we may use Eq.(l) plotting a(S,T)/ \ S — S c | p 
vs. T I | S — S c \ zu (not shown) with the three parame- 
ters S c , /i = v and z. The best data collapse is found for 
fi = 1.0 ± 0.1 and z — 2.94 ± 0.3, in very good agreement 
with the values obtained from Fig. 4. Additionally, we 
note the broad consistence with the earlier concentration 
tuning data where /i = 1.3 and z — 2.4 was inferred [p"4|| . 
We estimate the error of our combined analysis of the 
present stress-tuned data to 10% for /i and z. The criti- 
cal stress is determined with a relative accuracy to better 
than 0.1 kbar. It is important to note that either er(0) 
scaling (Fig. 2b) or dynamic scaling (Fig. 4) when taken 
by itself, may lead to a rather large error in /i and/or z, 
just because of the ambiguity of determining the critical 
region. However, the consistent determination of expo- 
nents from the combined scaling lends confidence to the 
values reported here. 

The above procedure to determine the conductivity 
at the critical stress is necessary because a c does not 
obey a simple power-law T dependence over the whole 
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FIG. 4. Scaling plot of a/a c vs. | S-S c \ / S c T y for sample 
1, with S c = 1.75 kbar and y = 0.34. 



T range. Above 100 mK the correction term dT w (with 
d < 0) comes into play. This is at variance with Si:B 
where a c ~ T 1 / 2 was observed in the whole range from 
60 to 800 mK (T|. On the other hand, a T 1 / 3 depen- 
dence of a in the vicinity of the MIT has been reported 
for transmutation-doped Ge:Ga over a large T range [ p^[ . 
Certainly, the finite-T behavior near the quantum criti- 
cal point needs closer theoretical scrutiny, in particular 
since the dynamic scaling is observed up to 800 mK when 
the correction term to cr c (T) is included. We remark that 
a simple algebraic T dependence a c = aT x which yields 
good dynamic scaling for Si:B p5[ , clearly leads to less 
satisfactory scaling in Si:P for any choice of x. 

In conclusion, we have demonstrated dynamic scal- 
ing of stress-tuned Si:P at the metal-insulator transition. 
The conductivity exponent fi « 1 is close to the expo- 
nents derived earlier from concentration tuning. How- 
ever, upon application of stress, the critical range is nar- 
rowed to conductivities below 12 _1 cm" 1 . Therefore, 
it is the absence of appreciable rounding effects in our 
samples close to the MIT that allows us to determine 
fi w 1 reliably, thus resolving the conductivity exponent 
puzzle. The temperature dependence of the conductiv- 
ity starting from the same cr(0) value for T — is dis- 
tinctly different for samples under zero stress and under 
stress. It is predicted that in the region between 15 and 
40f2 -1 cm -1 initially insulating stress-tuned samples will 
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exhibit a negative slope of cr(T), while samples under 
zero stress in this range are known to have a positive 
<y{T). In view of these differences away from the quan- 
tum critical point, the similarity of asymptotic dynamic 
scaling behavior is particularly noteworthy. A more de- 
tailed theoretical treatment which may eventually also 
account for the effective exponent fx sa 0.5 for samples 
above the crossover conductivity a cr is highly desirable. 
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